Fast Reconnection in High-Lundquist-Number Plasmas Due to the Plasmoid 

Instability 
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Thin current sheets in systems of large size that exceed a critical value of the Lundquist number 
are unstable to a super- Alfvenic tearing instability, referred to hereafter as the plasmoid instability. 
The scaling of the growth rate of the most rapidly growing plasmoid instability with respect to the 
Lundquist number is shown to follow from the classical dispersion relation for tearing modes. As a 
result of this instability, the system realizes a nonlinear reconnection rate that appears to be weakly 
dependent on the Lundquist number, and larger than the Sweet-Parker rate by nearly an order of 
magnitude (for the range of Lundquist numbers considered). This regime of fast reconnection is 
realizable in a dynamic and highly unstable thin current sheet, without requiring the current sheet 
to be turbulent. 
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The problem of fast reconnection in high-Lundquist- 
number plasmas has attracted a great deal of attention 
since the inception of the classical Sweet-Parker P, Q| 
and Petschek [3| models. In the Sweet-Parker model, the 
reconnection layer has the structure of Y-points, with a 
length of the order of the system size L, and a width 
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given hy Ssp = L/Sj^ , where Sl is the Lundquist num- 
ber based on the system size. The steady-state recon- 
nection time scale tsp is a strong function of Sl, and 
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increases as S"/ . For weakly collisional systems such as 
the solar corona, the Lundquist number is typically very 
large (~ 10^^ — 10^'') and hence, the time scale tsp is 
of the order of years, which is much too long to account 
for fast events such as solar flares. A similar difficulty 
appears when the Sweet-Parker model (equivalently, the 
Kadomtsev model is applied to the problem of saw- 
tooth crashes in high-temperature tokamaks, where the 
predictions of the model for the sawtooth crash time is 
significantly larger than that observed. These inadequa- 
cies of the Sweet-Parker model has led to the point of 
view that within the framework of the resistive MHD 
model, it is not possible to realize fast reconnection for 
high- 5 plasmas. 

In this Letter, we demonstrate that this widely ac- 
cepted point of view is not correct for large systems, 
characterized by high Lundquist number {Sl > 3 x 10^). 
By large systems, we mean systems for which the aspect 
ratio of the thin current sheet, Ssp/L, is much smaller 
than unity. In such systems, beyond a critical value 
of Sl, there appears to be three qualitatively distinct 
phases during the nonhnear evolution of reconnection. 
In the early nonlinear regime, the system evolves to a 
quasi-steady but transient state with the characteristics 
of Sweet-Parker reconnection — an extended thin current 
sheet in which reconnection occurs on the time-scale tsp- 
This quasi-steady regime is short-lived, and is followed 
by a second phase in which a rapid secondary instability 
(referred to hereafter as the plasmoid instability) occurs 
and produces plasmoids copiously. Although the pres- 
ence of this instability has been recognized for some time 
0, S 0, B B , it is only relatively recently that the pre- 



cise scaHng properties of the most rapidly growing plas- 
moid instability have been estabHshed[ll|, [IJ, llj. We 
demonstrate here that the key properties of this instabil- 
ity can be obtained as a simple extension of the classical 
dispersion relation of tearing modes applied to a 

Sweet-Parker current layer of width 6sp. The instability 
is followed by a third nonhnear phase in which the system 
exhibits rapid and impulsive reconnection mediated by a 
hierarchy of current sheets Q with widths that are much 
smaller than Ssp, with a Lundquist number dependency 
that decays faster than Sj^ When averaged over time, 
reconnection in the third nonlinear regime proceeds at a 
rate much faster than Sweet-Parker. Unlike the Sweet- 
Parker reconnection rate, which decays as Sj^ with in- 
creasing Sl, the reconnection rate in the third nonhnear 
phase appears to have a weak dependence on Sl- 

Linear Plasmoid Instability. To fix ideas, consider a 
Harris sheet of width a with the equilibrium magnetic 
field B — Bo tanh(z/a)x, which reverses sign at the z = 
surface. According to classical hnear tearing instability 
theory, the tearing mode growth rate for the " constant-')/'" 
and "non- constant-?/'" modes are as follows, respectively 
Q: 



(1) 

Here S = tr/ta, where ta — cl/Va — a{4:'Kpy^^ /Bq 
is the Alfven time, Tp = Aira^/rjc^ is the resistive diffu- 
sion time, p is the mass density, r] is the resistivity of the 
plasma, assumed to be a constant, and c is the speed of 
light. Note that the Lundquist number Sl, introduced 
earlier, is related to 5" through the relation, Sl — S/e, 
where e — a/L. According to the dispersion relation 
[11 the transition from "constant-i/i" to "non-constant-V'" 
modes occurs at kaS^^'^ ^ 1, with a peak growth rate 
that scales as ^maxTA S^^^"^. Since both branches of 
the tearing mode dispersion relation [T] yield instabilities 
that grow at a rate proportional to S raised to a frac- 
tional negative exponent, it may seem surprising upon 
first glance that there exists a super- Alfvenic tearing in- 
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stability that actually grows at a rate proportional to Sl 
raised to a fractional positive exponent. The key point is 
that unlike the Harris sheet, which is characterized by an 
equilibrium width a that is independent of the resistiv- 
ity, the Sweet-Parker current sheet is characterized by a 
width 6sp that has a strong dependence on the resistiv- 
ity. To connect with the results obtained in we in- 
troduce the dimensionless wave number k — KL = ka/e, 
and the dimensionless growth rate jl = JTA/e. In terms 
of dimensionless parameters k, e, and Sl, the maximum 
growth rate is given by jL.max ^ e^^^^K^^'^S^^^^ , and 
it occurs at the wave number Kmax ^ e^^^^S*^^^^. For 
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a Sweet-Parker current sheet, which obeys e ^ Sj^ , 

we then obtain Kmax ~ S^'^^i and the maximum growth 

rate ^L.max ^ S\l'^- Note that the number of plasmoids 
generated in the Hnear regime of this instability scales as 
i^max- These scaling relations are completely consistent 
with the results obtained for the fastest growing instabil- 
ity in 

The rapidity of this instability is not limited to Sweet- 
Parker current sheets, for we can easily extend the results 
obtained above to current sheets that may be narrower 
than Sweet-Parker current sheets (for the same value of 
Sl)- If we write e ~ SJ^" , where a > 0, we obtain 

^(5a-l)/4 jL,max ^ 5'(3a-l)/2^ folloWS then that 

in sheets for which a > 1/2, the plasmoid instability 
will grow even more rapidly than in Sweet-Parker current 
sheets of the same Lundquist number (for which a — 
1/2). 

Nonlinear Evolution of Extended Current Sheets. Our 
nonlinear study is carried in two-dimensional (2D), dou- 
bly periodic geometry using fully compressible MHD 
equations (in dimensionless form). The magnetic field 
is represented as B = VV'(a;, z) x y where V is a flux 
function and y is the ignorable coordinate, v = {vx,Vz) 
is the fluid velocity, and p is the scalar pressure, assumed 
to obey an isothermal equation of state, p = 2pT, where 
p is the density, and T is a constant temperature. The 
relevant dynamical equations are 

9tP + V-(pv) =0, (2) 
dtp + V ■ (pvv) ^ -Vp - VV'V^V, (3) 

dtip + v - ViIj = ryV^V- (4) 

Our numerical algorithm approximates spatial 

derivatives by flnite differences with a flve-point sten- 
cil in each direction, and time derivatives by a second- 
order accurate trapezoidal leapfrog method. We use 
variable grids in order to resolve the sharp spatial gra- 
dients in the reconnection layer. The initial condition 
consists of four flux tubes, which exhibit the coales- 
cence instability. The initial equilibrium state is given 
by ^p = Aq sin (27ra;/L) sin (27rz/L). In the results that 
follow, we take L = \/2 and Aq = L/2it, which produces 
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Figure 1: Time-sequence of the nonlinear evolution of the 
current density Jy of a Sweet-Parker current sheet in a large 
system of Lundquist number Sl = 6.28 x 10^ . The black lines 
represent surfaces of constant -if!. 



a peak magnetic fleld of unit magnitude in equilibrium. 
The initial density is taken to be approximately unity and 
the pressure (calculated from the isothermal equation of 
state with T — i throughout all simulations) is approxi- 
mately equal to 6. Both density and pressure profiles are 
assumed to have a weak spatial dependence in order to 
satisfy the magnetostatic equilibrium condition. 

The fiux tubes in equilibrium are unstable to the co- 
alescence instability, and produce an extended current 
sheet nonlinearly. Figure [1] shows a time-sequence of the 
current density Jy{x,z) for Sl = 6.28 x 10^, which ex- 
hibits the three stages of the nonlinear dynamics, dis- 
cussed above. (Here time is measured in units of tal — 
L/Va-) For convenience of visualization, we have rotated 
the current sheet by 45 degrees to make it horizontal. 
Frame (a) represents the first phase, the so-called Sweet- 
Parker phase of the instability, which is characterized 
by the formation of extended, system-size current sheets 
with the structure of Y-points. Subsequently, this cur- 
rent sheet becomes spontaneously unstable to the plas- 
moid instability, producing the second phase of nonlinear 
evolution (frame (b)). This second phase is followed by 
a third phase in which reconnection is typically impul- 
sive (or bursty), significantly faster than Sweet-Parker, 
and associated with the spontaneous development and 
ejection of plasmoids in the outfiow direction(s) (frame 
(c)). As mentioned above, the third phase has multiple 
stages. As the secondary instability develops, the initial 
Sweet-Parker current sheet breaks up, producing islands 
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Figure 2: Time-evolution of the scaled reconnection rate for 
three values of the Lundquist number. The black curve cor- 
responds to the case shown in Figure [T] 
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Figure 3: Averaged reconnection rate as a function of Sl Be- 
low Sl — ^ X 10'', the reconnection rate is Sweet-Parker and 

— 1/2 

decays as 5*^^ , shown by the dotted line. Above this thresh- 
old, the reconnection rate becomes weakly dependent on Sl. 



that separate shorter, narrower and more intense current 
sheets of width S < Ssp. These intermediate current 
sheets become unstable to secondary tearing instabihties 
of higher growth rate than those of the initial Sweet- 
Parker sheet, and lead to a regime of faster reconnection 
(frame (d)). 

In Figure [21 we show the time-history of the (normal- 
ized) reconnection rate for three simulations, including 
the high-^L simulation (black plot) shown in Figure [TJ 
The reconnection layer is along the x direction. We use a 
uniform grid along x. The grid along z is highly packed 
around z = 0. Table U summarizes the number of grid 
points as well as (the smallest) grid size in each direc- 
tion. (The Sweet-Parker thickness is included in for ref- 
erence.) The reconnection rate is calculated by taking 
the time derivative of the maximum stream function in 
the reconnection layer, normalized to the product of the 
upstream magnetic field and the Alfven speed, and mul- 
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tipHed by S^' . In Figure El the black curve, which corre- 
sponds to Sl = 6.28 x 10^, has been smoothed over OAta 



to eliminate rapid variations in a very dynamic current 
sheet. The linear instability of the initial state produces 
a Sweet-Parker current sheet that persists until about 
t ~ 3. After this point in time, the secondary instabil- 
ity develops, breaking up the long Sweet-Parker current 
sheet. (While the qualitative features of the instability 
are those of the plasmoid instability discussed above, it 
is difficult to test the precise scaling properties of the 
linear instability during a nonlinear evolution.) There is 
a rapid enhancement in the reconnection rate that re- 
mains impulsive and ffuctuates widely even as the rate 
tends to settles down to a plateau, until about t ~ 9. 
At this stage of the third nonlinear phase, some of the 
small islands produced by the secondary instability coa- 
lesce to form larger islands that are convected towards the 
boundaries. (If the islands grow to large size but are con- 
strained to stay fixed at the center of the computational 
domain by reason of symmetries imposed in the simula- 
tions, the third nonlinear phase may be short-lived, and 
the reconnection rate may fall rapidly.) At about this 
point in time, the extended current sheet shows yet an- 
other burst of secondary tearing activity producing mul- 
tiple plasmoids, and a consequent enhancement in the 
reconnection rate, which at about t ~ 12 attains an or- 
der of magnitude higher than the Sweet-Parker rate at 
this value of Sl- Due to insufficient spatial resolution, 
caused by the slow drift of the current sheet away from 
the region where the grid points along z are clustered, we 
are not able to carry these simulations forward longer in 
time. 

The plasmoid instability of Sweet-Parker sheets occurs 
after Sl exceeds a critical value, determined numerically 
to be approximately 3 x 10'' in the present study. Like 
the black curve, the blue dashed curve in Figure [21 corre- 
sponds to another value of Sl {— 2.51 x 10^) above the 
threshold and shows generically similar behavior, while 
the red dashed curve corresponds to a value (= 3.14x 10^) 
at about the threshold. We have identified the plateau 
regions in these curves by inspection, and stated the cor- 
responding time-intervals in Table [H There is clearly an 
element of arbitrariness in this identification, because the 
underlying dynamics is not quasi-steady. In fact, the dy- 
namics is generically impulsive, strongly time-dependent, 
involving the growth and ejection of multiple plasmoids 
in the presence, and complicated by the sloshing of co- 
alescing islands as well as the tilt and slow drift of the 
current sheet away from the central region. The slosh- 
ing is most clearly seen in the oscillations of the red 
curve, where the first peak corresponds to the highest 
reconnection rate, after which the system settles down to 
a Sweet-Parker-like plateau regime. A rough definition 
of the plateau region is when the reconnection rate ap- 
pears to fiuctuate about a relatively steady value. If the 
Lundquist number of the simulation is above the thresh- 
old for the plasmoid instability, the plateau region is typ- 
ically characterized by the presence of multiple islands. 

Figure [3| shows a plot of the averaged reconnection rate 
(normalized) as a function of Sl- The averaged recon- 
nection rate is calculated by averaging the smoothed re- 
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connection rate curve over the plateau region. The error 
bar is calculated from the standard deviation of the data 
points in the plateau region. Our simulations suggest 
that below a threshold, S'l = 3 x 10'*, the reconnection 
rate exhibits Sweet-Parker scaling. Over this threshold, 
the averaged reconnection rate exhibits a much weaker 
dependency on Sl- At high values of Sl, the average re- 
connection rate exceeds the Sweet-Parker rate by nearly 
an order of magnitude. We conjecture that as the initial 
Sweet-Parker current sheet is broken up into sufHciently 
short segments, and the effective Lundquist number of 
these segments falls below the threshold of instability, 
the process will attain a quasi-steady state, provided the 
reconnection dynamics is not quenched by the erosion 
of the flux content of the system due to the cumulative 
effects of resistive diffusion. 

These results bear qualitative resemblance and are 
complementary to those reported in a concurrent study 
[lB|, where 2D, fully electromagnetic particle- in-cell sim- 
ulations that include a collision operator show a regime 
of fast reconnection brought about by the onset of sec- 
ondary tearing instabilities. Since the physical models 
underlying the simulations presented here and in are 
quite different, it is not surprising that there are some im- 
portant differences between the two studies. The initial 
Sweet-parker current sheet realized in shows clear 
evidence of the plasmoid instability, and the width of 
subsequent intermediate current sheet segments becomes 
much narrower as time progresses, as in the present study. 
Whereas this narrowing is controlled entirely by the value 
of the Lundquist number and can continue without the 
intervention of kinetic effects in the present study, ki- 
netic effects inevitably intervene in when the current 
sheet becomes of the order of ion skin depth. For sys- 
tems of moderate size (of the order of 200di or less) , when 
S < di, the simulations in [lB| show the onset of fast re- 
connection, consistent with the predictions of earlier Hall 
MHD models [ITi, [13, 120] , and verified by experiment 
[2l] |. However, for systems of larger size, the transition to 
fast reconnection appears to occur at values of Lundquist 
number that are somewhat smaller than the Hall MHD 
prediction [TtI . [T§ . [l9l . [20| due to the intervention of the 
plasmoid instability. 

In summary, we have demonstrated that within the 
framework of the resistive MHD model, thin current 
sheets in systems of large size that exceed a critical value 
of the Lundquist number are unstable to a super- Alfvenic 
tearing instability, referred to here as the plasmoid insta- 
bility. If the aspect ratio of the current sheet is written 
e ~ S^"', where a > 0, the fastest growing instability is 

characterized by a wave number Kmax ^ S^^" ^■'^^ and 

a growth rate jL.max sf" In the special case 

of a Sweet-Parker current sheet for which a = 1/2, we 
obtain the instability discussed in [ll|. As a result of this 
class of instabilities, the system realizes a nonlinear re- 
connection rate that appears to be weakly dependent on 
the Lundquist number, and larger than the Sweet-Parker 
rate by an order of magnitude (for the range of Lundquist 



Table I: Details of the spatial resolution and the plateau 
regime for the runs shown in Figure [2] 
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Az 








Plateau Regime 


3.14e4 
2.51e5 
6.28e5 


6.7e-4 
5e-4 
3.3e-4 


5e-4 
9.4e-5 
3.9e-5 


1501 
2001 
3001 


501 
1001 
1001 


5.6e-3 
2e-3 
1.26e-3 


2.9-8.9 
6.9-13.0 
6.9-11.6 



numbers considered in the present study). We note that 
this weak dependency on the Lundquist number is real- 
ized in a thin current sheet that is violently unstable and 
dynamic, but not turbulent. In this respect, our results 
are different from those obtained from turbulent recon- 
nection studies in 2D [I^, [H, as well as 3D [H,!!!!, 
which also report reconnection rates that depend weakly 
on the resistivity. 
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